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XMM-NEWTON OBSERVATIONS OF IGR J0029 1+5934: SIGNS OF A THERMAL SPECTRAL 

COMPONENT DURING QUIESCENCE 
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ABSTRACT 

We present X-ray observations of the transient accretion-powered millisecond pulsar IGR J0029 1+5934 during 
quiescence. IGR J00291+5934 is the first source among accretion powered millisecond pulsars to show signs 
of a thermal component in its quiescent spectrum. Fitting this component with a neutron star atmosphere or a 
black body model we obtain soft temperatures 64 eV and ^110 eV, respectively). As in other sources of 
this class a hard spectral component is also present, comprising more than 60% of the unabsorbed 0.5-10 keV 
flux. Interpreting the soft component as cooling emission from the neutron star, we can conclude that the compact 
object can be spun up to milliscond periods by accreting only ^ 0.2 Mq. 

Subject headings: accretion, accretion disks — star: individual (IGR J00291+5934) — stars: neutron 



1. INTRODUCTION 

IGR J0029 1+5934 (IGR J0029 in the following) is one out 
of eight low mass neutron star transient showing coherent pul- 
sations during outbursts (see Wijnands 2005 for a review). 
Distinctive properties of these faint transients are weaker out- 
burst peak luminosities (~ (1 — 5) x 10^^ erg s~^) with re- 
spect to classical neutron star transients (see Campana et al. 
1998 for a review), very low mass companions (mass functions 
< 2 X 10^^ Mq), short orbital periods (Poib ^ 4.3 hr), very 
faint quiescent X-ray luminosities ( ;S 5 x 10^^ erg s~^) and 
absence of a soft X-ray spectral component in quiescence (qui- 
escent spectrum consistent only of a hard power law compo- 
nent). 

IGR J0029 was discovered while in outburst on December 
2nd 2004 during a routine monitoring of the Galactic plane with 
the INTEGRAL satellite (Eckert et al. 2004). Follow-up RXTE 
observations revealed 598.88 Hz (1.67 ms) coherent pulsations 
(Markwardt, Swank & Strohmayer 2004), making it the fastest 
confirmed accretion-powered millisecond pulsar and a 2.46 hr 
orbital modulation (Markwardt et al. 2004). In quiescence IGR 
J0029 was detected by Chandra at the beginning of 2005 (and 
historically by ROSAT; Jonker et al. 2005). The source quies- 
cent 0.5-10 keV unabsorbed flux was ^ 8 x 10^^'* erg cm^^ 
s~^ with a variability among different observations by a factor 
of 2 (Jonker et al. 2005). Chandra detected IGR J0029 also 
in November 2005 at a level of 7 x lO"^** erg cm^^ s^^ 
(Torres et al. 2008). Finally a short observation was carried 
out with Chandra in September 2006 with the source at a level 
of - 1 X IQi^ erg cm'^ ^-i (jonker, Torres & Steeghs 2008). 
All these Chandra observations collected a relatively low num- 
ber of photons and a detailed spectral analysis could not carried 
out, being all spectra consistent with a single power-law (pho- 
ton index 2.4 ± 0.5, 90% confidence level). 

The non-detection of X-ray bursts makes the estimate of the 
distance to IGR J0029 difficult. Shaw et al. (2005) gave an up- 
per limit of 3.3 kpc considering the source position with respect 
to the Galactic Center. Torres et al. (2008) estimated a distance 
of 2 — 4 kpc based on the critical X — ^ray luminosity needed to 
ionize the accretion disc and produce the observed X-ray light 
curve during outburst. We will assume in the following a fidu- 



cial distance of 3 kpc. 

In this letter we report on the first study of IGR J0029 in 
quiescence with XMM-Newton. The larger throughput of the 
XMM-Newton instruments allowed us to study in detail the qui- 
escent spectrum. In Sect. 2 we discuss the XMM-Newton data 
analysis. Sect. 3 contains our discussion and conclusions. 

2. DATA ANALYSIS 

The XMM-Newton Observatory (Jansen et al. 2001) in- 
cludes three 1500 cm^ X-ray telescopes each with a European 
Photon Imaging Camera (EPIC, 0.1-15 keV) at the focus. Two 
of the EPIC imaging spectrometers use MOS CCDs (Turner et 
al. 2001) and one uses pn CCDs (Struder et al. 2001). 

XMM-Newton observed IGR J0029 on July 24 2007 for 27 
ks with the thin filter on all the EPIC instruments. Data were 
processed using SAS v.6.6.0. At the end of the observation a 
strong soft proton flares occurred limiting the usable observing 
time to 17 and 24 for the pn and the MOS detectors, respectively 
(filtering out background flares for total 0.2-15 keV rates less 
than 5 and 20 c s^^ for the MOS and pn cameras, respectively). 
We extracted the pn source spectrum from a circular region with 
radius 22" and the background spectrum from a larger region 
close to the source with radius 34" within the same CCD. For 
the MOS cameras we used the same source extraction region 
and a 84" background region. We obtained 466, 142 and 146 
counts for pn, MOSl and MOS2 cameras in the 0.2-10 keV 
range. The background contribution is at a level of 50% and 
30% for the pn and MOS cameras, respectively. 

We considered for timing analysis only the pn light curve in 
consideration of the higher signal to noise ratio. We subtracted 
the background light curve and we did not find any evidence 
for variability down to the 500 s, i.e. the shortest timescale 
that could be investigated (the fit with a constant provides a re- 
duced = 0.8 with 43 degrees of freedom, dof). This is at 
variance with what observed in the optical band, where strong 
flares (~ 1 mag) were reported on timescale of >, 500 s (Jonker 
et al. 2008). 

The spectral analysis was carried out by using the data 



from the three EPIC cameras together. Spectra were binned 
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to, at least, 20 counts per spectral bin resulting in 37 bins. 
A single power law fit provides a good description of the 
data (reduced = 0.9, 34 dof) but the inferred column 
density (using TBABS model in XSPEC12.4.0ai) is Nh = 
1.61 X 10^^ cm ^ (90% confidence level), which is much 
smaller than the column density observed in outburst of Nh = 
4.3to5 X 10^^ obtained fitting together Chandra HETG 
and RXTE/PCA data (Paizis et al. 2005). This value is also 
consistent with the Galactic value of Nh = 4.6 x 10^^ cm~^. 
Here we adopt this value (see Table 2). 

Fitting the data with the column density fixed at the above 
value provides relatively poor fits (reduced = 1-35 in the 
case of a power law). Therefore we tried double component 
models. Usually spectra of neutron star transients in quiescence 
are fit with a two components model: a neutron star atmosphere 
emission model (or a black body) plus a hard power law com- 
ponent. These models provide a better description of the data. 
An F-test gives a 3.4 fj probability for the significance in the fit 
improvement due to the extra model component (either black 
body or neutron star atmosphere, see Table 2). This indicates 
that, if Nh is held fixed at the value measured in outburst, the 
addition of a soft spectral component is significant. Estimates 
of the column density during quiescence of transient accreting 
milUsecond X-ray pulsars are available only for SAX J 1808. 4- 
3658. For this source the quiescent value of the column density 
is consistent with the one observed in outburst (Heinke et al. 
2007; Campana, Stella & Kennea 2008). 

3. DISCUSSION AND CONCLUSIONS 

In recent years the quiescent properties of a handful neutron 
star transients have been studied in some detail thanks to dedi- 
cated observations mainly carried out with XMM-Newton and 
Chandra. 'Classical' (i.e. not showing pulsations during out- 
bursts) neutron star transients have quiescent (0.5-10 keV) lu- 
minosities within a narrow range of 10^^ — 10^^ erg and 
their spectra display a soft component consistent with (cool- 
ing) emission from the entire neutron star and in some cases a 
hard (power-law) component contributing 10 — 50% of the flux 
(e.g. Campana et al. 1998). An increasing number of candidate 
neutron star systems has also been discovered in globular clus- 
ters (e.g. Pooley et al. 2003; Heinke et al. 2003). These are 
characterized by a soft spectrum consistent with emission from 
the entire neutron star surface, but do not display any evidence 
for a power law component (upper Hmits on the flux of such a 
component are very tight in some cases). No outbursts have yet 
been detected from these sources. At the other extreme, there 
exist systems like PSR J0024-7204W in the globular cluster 
47 Tuc, showing emission consisting of a dominant power law 
component, which is eclipsed for a portion of the orbit, and a 
soft thermal component, which appears to be persistent (Bog- 
danov, Grindlay & van den Berg 2005). 

Transient accreting millisecond X-ray pulsars have low qui- 
escent (0.5-10 keV) luminosities, below (and in some cases 
well below) 10^^ erg s~^. In these systems the quiescent spec- 
trum is made just by a power law component with no signs of a 
soft component. The best known system, SAX J1808.4-3658, 
shows a pure power law spectrum, with an upper limit on the lu- 
minosity of the soft component of ^ 23% (Heinke et al. 2007). 

Our XMM-Newton data on IGR J0029 in quiescence pro- 
vides for the first time evidence for the presence of a soft 



component in the quiescent spectrum of a transient accreting 
millisecond X-ray pulsar. The temperature of this component 
(fc T ~ 70 eV, when modeUed with a neutron star atmosphere) 
is soft, but in line with other classical neutron star transients 
(e.g. Heinke et al. 2007). The bolometric neutron star atmo- 
sphere luminosity is 1.4 x 10"^^ erg s~^. This luminosity is 
usuaUy ascribed to the heat generated deep in the crust during 
outburst episodes (i.e. accretion onto the neutron star surface) 
which is then radiated away through the atmosphere while in 
quiescence producing the quiescent X-ray emission (Brown et 
al. 1998; Haensel & Zdunik 1990). The efficiency of die re- 
radiation process depends on which fraction of heat escapes the 
neutron star as neutrinos rather than photons and in turn on the 
physical conditions (i.e., composition, density, and pressure) of 
the neutron star interior 

The accretion history of IGR J0029 is relatively well known. 
IGR J0029 was discovered during an outburst in Dec. 2004. 
The bolometric fluence of this outburst has been estimated in 
1.8 X IQ-^ erg cm^^ ^sj„g rxTE data (Galloway et al. 2005) 
or 1.37 X 10"^ erg cm'^ using INTEGRAL data (Falanga et 
al. 2005). Thanks to file RXTE-ASM two previous instances 
of activity, on 1998 November 26-28 and 2001 September 11- 
21 have been revealed (Remillard 2004). These results indicate 
that the transient appears regularly, roughly every 3 yr. Indeed a 
new outburst started on Aug 13 2008 (Chakrabarty et al. 2008), 
providing further evidence that the outbursts recur fairly regu- 
larly. Assuming that all the outbursts have comparable fluence 
(as testified by similar outbursts peak in the RXTE- ASM), one 
can estimate a long-term J'acc = 1-4 x 10~^^ erg cm~^. 
Based on deep crustal heating theory, the quiescent flux is re- 
lated to the long term time-averaged flux Fq = i^acc/135 = 
1.0 X 10"^^ erg cm~^ (Brown et al. 1998). This is in line 
with what observed. The goodness of this match lends support 
to the idea that the soft spectral component derives from cool- 
ing of the neutron star surface and not from accretion (since the 
spectrum at these low luminosities is expected to be soft, e.g. 
Zampieri et al. 1995). 

Having estimated the quiescent cooUng luminosity and the 
time-averaged flux (or mass inflow rate) we can put constraints 
on models for the neutron star interior (Yakovlev & Pethick 
2004; Levenfish & Haensel 2007). As can be seen from Fig. 
2, IGR J0029 lies on the standard cooling curve for a low mass 
neutron star. This is at variance with, e.g. SAX J1808.4-3658 
and Cen X-4, for which low cooUng luminosities have been 
inferred, hinting for enhanced neutrino emission produced in 
the high-density core of a relatively high mass neutron star 
(M > 1.6 - 1.7 Mq, e.g. Colpi et al. 2001). 

Our observation indicates that neutron stars showing coher- 
ent pulsation during their outbursts do not necessarily lack a 
soft component in their quiescent spectra. If the soft spectral 
component is due to deep crust heating emission and the pres- 
ence or lack of this component is a proxy of its mass, as mas- 
sive neutron stars do not show it due to fast cooling. This would 
imply that light (M ^1.6 Mq) neutron stars too can show co- 
herent pulsations at the milhsecond level, having accreted only 
< 0.2 Mq (Lavagetto et al. 2004). 



SC thanks C. Heinke for providing the SM macro and data 
for producing Fig. 2. 
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Fig. 1 . — XMM-Newton spectrum of IGR J0029. In the upper panel data along with the best fit (absorbed) NSA plus power law spectrum are shown: pn data 
are marked with a open squares, MOSl data with open circles and MOS2 data with filled circles. In the bottom panel the contribution of each bin. Dashed and 
dotted lines refer to the neutron star atmosphere and power-law spectral components, respectively. 



Table 2 
Spectral fits. 



Model 


Nh 
(lO^i cni-2) 


Ph. Index 


Temperature 
(keV/log T) 


Xrcd 

(dof) [nhp°] 


Flux'' 
(Luminosity) 


Power law 
Black body 
Power law 
Black body 

USA" 
PL+BB'^ 
PL+NSA'^ 


1 6+*^'' 

^■^-0.8 

< 0.2 
4.6 
4.6 
4.6 
4.6 
4.6 


1 7+U.2 
^- ' -0.2 

9 c+O.S 
^■3-0.3 

1 e+0.3 
-■^•^-0.5 


"•3O-0.07 

^■^^-0.03 

t- 07+0.03 


0.89 (34) [66] 
1.55 (34) [2] 
1.35 (35) [8] 
2.86 (35) [0] 
3.46 (35) [0] 
0.93 (33) [59] 
0.98 (34) [49] 


7.2 (7.7) 

3.4 (3.6) 

8.5 (9.2) 
4.4 (4.7) 
7.1 (7.7) 

12.4(13.3) 
11.6(12.4) 



All quoted errors are at 90% confidence level (Ax = 2.71). 



"Null hypothesis probabihty. 

''Unabsorbed 0.5-10 keV flux in units of 10"^'* erg cm~^ s~^ and luminosity in units of lO^'^ erg s~^ at 3 kpc. 

'^NSA normalization fixed at 1.1 X 10^^ as for a source at 3 kpc. Neutron star mass, radius and magnetic field fixed to 1.4 Mq, 10 km and G, respectively. 
'^Black body radius R = y km at a distance of 3 kpc. The black body component comprises 36% of the 0.5-10 keV unabsorbed flux. 

"^The NSA component comprises 38% of the 0.5-10 keV unabsorbed flux. 



4 



— \ — I — I — I — \ — I — I 

Aql X-1 

NGC 6440 

4U 2129+47 

KS 1731-260 

4U 1608-522 

Terzan 5 

IM 1716-315 

RX J1709-26 

MXB 1659-29 
XTE J2123-058 
SAX 1810.8-2609 
IH 1905+000 
2S 1803-245 



~| — \ — I — I — I — I — I — I — r 

A XTE J0929-314 

B SAX J1808. 4-3658 

C XTE J1807-294 

D XTE J1751-305 

E XTE J1814-338 

F IGR J00291+5934 



9 

10 
11 
12 
13 



Dxl.5 




-11 -10 
log M (Mg/yr) 



Fig. 2. — Cooling curves for various neutron star interior neutrino emission scenarios, compared with measurements (or 95% contidence upper limits) of the 
quiescent bolometric luminosity and time-averaged mass transfer rate for several neutron star transients (taken from Heinke et al. 2008). The big (green) dot 
indicate IGR J0029 values as estimated in this paper Uncertainties have been accounted for by using a distance range of 2-A kpc. The cooling curves are taken from 
Yakovlev & Pethick (2004); the dotted curve represents a low-mass neutron star, while the lower curves represent high-mass NSs with kaon or pion condensates or 
direct Urea (Durca) processes mediated by nucleons or hyperons. 
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